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Abstract. The pyrochlore antiferromagnet Tb2Ti207 has proven to be an 
enigma to experimentalists and theorists working on frustrated magnetic systems. 
The experimentally determined energy level structure suggests a local (111) Ising 

antiferromagnet at low temperatures, T ~ lOK. An appropriate model then 
predicts a long-range ordered Q = state below approximately 2K. However, 
muon spin resonance (fiSH) experiments reveal a paramagnetic structure down 
to tens of milli-Kelvin. The importance of fluctuations out of the ground 
state effective Ising doublet has been recently understood, for the measured 
paramagnetic correlations can not be described without including the higher 
crystal field states. However, these fluctuations treated within the random phase 
approximation (RPA) fail to account for the lack of ordering in this system below 
2K. In this work, we briefly review the experimental evidence for the collective 
paramagnetic state of Tb2Ti207. The basic theoretical picture for this system 
is discussed, where results from classical spin models are used to motivate the 
investigation of quantum effects to lowest order via the RPA. Avenues for future 
experimental and theoretical work on Tb2Ti207 are presented. 



1. Introduction 

Frustrated magnetism has been actively studied for decades, ^ |21 EI and interest 
continues to grow. In models like the antiferromagnetic (AFM) J1-J2 model on a 
bipartite lattice, frustration results from a competition between two different ground 
state spin structures, one favored by Ji and the other by J2. 0] The J1-J2 model 
has motivated many interesting theoretical studies, but the tuning of energy scales to 
induce frustration is not easily achieved experimentally. Another class of materials 
and models exhibit frustration at the level of nearest neighbor interactions via a 
competition imposed by the lattice geometry, creating geometric frustration. [2j 
Crystal geometries formed from unit triangles and AFM interactions are frustrated 
because the spins about a triangle can not arrange themselves so that all pairwise 
interactions are satisfied (i.e., — JiS^ • Sj for Ji < can not be minimized for each pair 
Hence, the triangular and kagome lattices are common frustrated geometries 
in 2d. In three dimensions, lattices of corner sharing triangles (garnet structures) 
or tetrahedra (pyrochlore and spinel structures) are frustrated with nearest neighbor 
AFM interactions. 
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The pyrochlore lattice has received increased attention recently because the 
5=1/2 AFM model is predicted to be a collective paramagnet, or spin hquid, 
in the ground state. [HI E] A 3(i spin liquid is unusual, [S| but the role of lattice 
dimensionality appears to be reduced in the face of large quantum fluctuations driven 
by frustration. [Hj The stabilization of a non-coUinear long-range ordered state is 
expected when long-range interactions are included, 101 or when constraints arising 
from single-ion anisotropy are imposed. 10^ Magnetic insulating pyrochlores, general 
formula A2B2O7, where A is a rare earth ion (Ho^+, Dy'^"'", Tb'^"'', Gd'^+), are a subset 
of frustrated magnets that are at the focus of much experimental and theoretical 
work. 012] This family of materials displays long-range order, JI] a novel "ice like" 
phase with residual entropy, jH^l^] and collective paramagnetic behavior at milli- 
Kelvin temperatures, jj^ The focus of our brief review is the candidate 2>d spin liquid 
Tb2Ti207. |14j The many conflicting experimental results and theoretical predictions 
of this material are discussed, and the recent progress toward a general picture for 
Tb2Ti2 07 is presented. 

2. Experimental Picture of Tb2Ti207 

The pyrochlore Tb2Ti2 07 has dominant AFM interactions as determined from dc 
magnetization measurements on a polycrystalline sample, 6cw ~ — 20K. ^3] Results 
from muon spin resonance (/xSR) clearly indicate dynamic moments down to 70mK 
despite the development of short-range correlations at approximately 50K, as seen in 
neutron scattering measurements. Suppression of an ordering temperature 

in this material by two orders of magnitude (i.e., |0cw|/7OmK = O(IO^)) can be 
understood in the context of a nearest neighbor AFM. However, Tb2Ti207 has a 
measured spin anisotropy gap of A w 18K, ^^E] attributable to crystal field effects, 
as well as reasonably strong dipolar interactions, EI that such a large reduction 
in Tc is quite puzzling. 

As a first step to unraveling the mystery of Tb2Ti207, one has to determine 
the single ion properties of Tb'^"'" in the dense material. From experiments and 
calculations, Gingras et al. find that the Tb'^^ moment is approximately 5/iB, a 
substantial reduction from the free ion result 9.6/iB. 16 In addition, the crystal field 
effects on Tb^+ (5* = 3, L = 3 so J = 6) yield an effective Ising ground state level 
structure with a quantization axis oriented along the local (111) direction, 16, 18 with 
a gap to the next energy levels, an excited doublet, on the order of A « 18K. This 
(111) anisotropy is significant since for AFM effective nearest neighbor interactions 
it greatly reduces the frustration from that of a Heisenberg like model. JLO, With 
a better grip on the dipolar strength and the crystal field levels, a revised Curie- 
Weiss temperature of 6'cw ~ — 14K for the dipoles and exchange is calculated and an 
ordering temperature of IK deduced. jl6lll7j Also, from theoretical work on a dipolar 
spin ice model ((111) Ising) one is able to show that for the parameters of Tb2Ti207, 
a long-range ordered state with zero net moment about each tetrahedron (an all-in 
all-out state) is expected at Tc ~ I.IK. ^7] That experiments prove these predictions 
wrong is possibly due to the effects of long-range exchange beyond nearest-neighbor 
or quantum fluctuations. 

With the single ion properties suggesting a (111) Ising dipolar model for Tb2Ti207 
at T < lOK, but the absence of any agreement between this model and experiments 
on an ordering temperature, a re-evaluation of the model is necessary. Hence, the 
focus of the discussion shifts to the paramagnetic (PM) regime and the short-range 
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correlations. Early inelastic neutron data on a single crystal found evidence for the 
partial softening of the lowest energy modes about [002] and [220] . ^1 Elastic 
scattering measurements on a single crystal of Tb2Ti2 07 clearly display diffusive 
scattering about these same regions in the (hhl) plane with the strongest correlations 
centered about [002]. ^11201 The qualitative features of the structure factor (^(q)) in 
the PM phase can be roughly fitted with a near neighbor AFM Heisenberg model, fTHj 
suggesting that the strict Ising anisotropy might not hold for temperatures defining 
the PM regime of this material, i.e., T ^ 0.5K. Note, the PM S'(q) is measured at 9K 
in Gardner et al. but at 0.4K in Yasui et al. pDj A relaxation of the (111) spin 
anisotropy in a long-range dipolar AFM pyrochlore model does yield good agreement 
with the PM S'(q) in this material for temperatures T - \ecw \ ^ lOK. [HIES 

Recent data for Tb2Ti207 in the low temperature regime suggests the existence 
of a spin glass like state. Hysteresis is observed in the field cooled and zero field cooled 
neutron scattering data below 1.5K. [521 Thermodynamic measurements of the static 
susceptibility also indicate spin glass physics, albeit at a much lower temperature, 
70mK. I2SI- 

Hence, at present the general consensus is that Tb2Ti207 remains a collective 
paramagnet down to T = 50-70 mK at ambient pressures. Under high pressure, 2 — 8 
GPa, the material orders at Tc « 2K, |22 but the spin structure is not the Q = 
all-in all-out structure predicted from a (111) Ising model, ^| 122] or the "nearby" 
long-range ordered spin ice state. j2^ An applied field can also induce an ordered 
state in this material. Field driven ordered states is a current topic in the spin 
ice material Dy2Ti207, 27 but unlike the spin ice materials moderate field strengths 
along high symmetry directions indicate weaker spin anisotropy for Tb2Ti2 07. |20| 
We note that the gap between the ground state doublet and the first excited states 
in Tb2Ti207 is an order of magnitude smaller than the gap measured in the spin ice 
materials. ^HI 

3. Theoretical Interpretation of Tb2Ti207 

The Heisenberg AFM on the pyrochlore lattice is defined by the simple Hamiltonian 
H = —JiJ2{{i a) {j b)) ' ^5 ' where Ji < 0, (i,a) describes the location of a spin 
by its fee lattice position Ri and tetrahedral basis point r°, and interactions are 
among nearest neighbors only. The Hamiltonian H represents a highly frustrated 
model that produces a spin hquid ground state for quantum spins (Sf = 1/2) [HI [7] 
or classical spins (S^ = xS^''' + ySf'^ + zSi'""). EH |^ EOl ED In the classical model, 
an ordered state can be induced by perturbations arising from exchange interactions 
beyond nearest neighbor jSHI or by constraints imposed by spin isotropy. jTH] If the 
single-ion anisotropy is along the local (111) direction (i.e., = z°-af where i° is 
the local (111) direction and af — ±1), then it is the nearest neighbor ferromagnetic 
pyrochlore (nearest neighbor spin ice) that is frustrated. [21 El E] For rare earth 
magnets like the spin ice materials and Tb2Ti207, long-range dipolar interactions are 
a significant contribution to the Hamiltonian. Therefore, the basic model for rare 
earth pyrochlore insulators can be written. 




3(Sf ■R^)(S;.R°f) 
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where D^^ ~ ^^nn the nearest neighbor dipole strength, D — fioji^ /4:TTRf-^^, ^ is 
the moment on the rare earth ion, is the nearest neighbor distance, and R"^'' is 
the vector separation between spins and Sj. If the spins in Eq.QJare constrained to 
their local (111) Ising axes, then Eq.^ represents the dipolar (111) Ising model (DIM) 
that successfully describes the specific heat JT] and paramagnetic correlations 32 of 
spin ice. When the DIM is simulated by a Monte Carlo algorithm employing nonlocal 
loop moves, a first order transition to a spin ice ground state is found, j25j with an 
ordering wave vector and spin structure in agreement with mean-field calculations. j21| 

The phase diagram of the DIM on the pyrochlorc lattice is generated by defining 
the sign and magnitude of Ji and the strength of the dipole-dipole interaction D. |17j 
The DIM exhibits two phases depending on the ratio Ji/D: For Ji/D ^ —4.525 a spin 
ice manifold results (macroscopic degeneracy) , where two spins point in and two spins 
point out along their local (111) axes for each unit tetrahedron. For Ji/D < —4.525, a 
non-collinear long-range ordered state results, where all spins about a unit tetrahedron 
point either in or out along their local (111) axes, this Q = structure is replicated 
over the lattice. We note that going from the PM to the Q = state occurs via 
a second order phase transition, but the path from the PM to spin ice manifold 
represents a cross over to a dynamically frozen state without long-range order. In the 
DIM, the temperature of the phase transition or cross over is signalled by the peak in 
the calculated specific heat, whose position depends on the value of Ji/D. |E| Hence, 
the important energy scale in this model is the effective nearest neighbor exchange, 
Jeff = Ji/3 + 515/3, where 1/3 and 5/3 are geometric factors arising from the inner 
product of the local (111) spin quantization axes. With Ji/D « —5.5 for Tb2Ti207, 
the prediction is a Q = state at approximately IK. |lf)lll7j The subtle point about 
Tb2Ti207 in the DIM is that it sits near the boundary between the two phases, 
making it susceptible to fine tuning of the model parameters (exchange or dipolar 
energy scales and single ion crystal field state wave functions). It is conceivable that 
Tb2Ti207 might sit close to a line of disorder in the phase diagram where fluctuations 
dominate and a tendency to order is suppressed. The simulation of the DIM near the 
phase boundary presents some numerical challenges and is not well understood. j25| 

As mentioned in Section|21 Tb2Ti207 does not order at T w IK, and suggestions 
from the PM S'(q) results ^11201 that the (111) Ising constraint should be relaxed 
in this case provides a hint of where to venture next. In considering the symmetry of 
the DIM, one can show that no strict (111) Ising model on the pyrochlorc lattice can 
reproduce the short-range correlations observed in experiments, but a model in which 
the rotational symmetry is partially restored is found to admit a solution compatible 
with experiments. ,21, Mean-field and Monte Carlo calculations make these claims 
more quantitative. Hence, in the PM regime the physics of Tb2Ti207 can not be 
ascribed to the same (111) Ising interpretation of the moments that works so well for 
the spin ice materials. 

A pitfall with a classical Heisenberg description for Tb2Ti207 is that complete 
restoration of spin isotropy does not agree with experiments and calculations of the 
single ion properties of the Th^^ . [TB1[TH] The g-tensor of Th^+ in Tb2Ti207 exhibits 
(111) like anisotropy below lOK and a ground state doublet structure that is separated 
from a first excited doublet by about 18K, which is on the order of the Curie- Weiss 
temperature for exchange and dipolar interactions (^cw — — 14K). jl()| For the spin 
orbit coupled Tb'^^ moments (J — 6), the angular momentum for the two lowest 
lying doublets is predominately \Mj = ±4) for ground state and \Mj = ±5) for the 
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first excited state. jTC] Hence, transverse fluctuations between these two lowest lying 
doublets are allowed and non-negligible at the temperatures of the S{c\) experiments, 

5-ioK. f^iEa 

Applying a random phase approximation (RPA) to a model with exchange and 
long-range dipoles, the fluctuations between the lowest lying crystal field states can 
be studied, i.e., the (111) Ising limit is recovered when transitions between states 
IV'o) (ground state) and l^/^i) (first excited state) have a vanishing matrix element, 
(V'll^^lVo) = 0. Using |i/.o) = \Mj = ±4) and |Vi) = \Mj = ±5 >, the 
RPA 122 gives excellent semi-quantitative agreement with the PM ^(q) and the energy 
dispersion between the low lying magnetic states of the neutron experiments. |19| 
However, the RPA still predicts an ordered state at Tc « 2K. Hence, the RPA 
demonstrates the sensitivity of the PM correlations observed in neutron scattering 
to quantum fluctuations in this unusual magnetic system (|'0o) = \Mj = ±4) and 
I'i/'i) = \Mj ~ ±5 >), but it still gives a Tc of approximately 2K to the above all-in 
all-out Q = state. 

4. Avenues for Future Studies and Conclusions 

In order to understand the physics of Tb2Ti207 below 2K, we need to go beyond the 
first order approximation of the RPA and consider quantum effects more rigorously. 
A reasonable first step would be to employ a more detailed description of the crystal 
field states of the Tb3+ ion, and then form an effective low energy Hamiltonian 
from the more complicated model following a procedure similar to that used in the 
derivation of the t-J or Heisenberg models from the Hubbard model. Experimentally, 
low temperature neutron scattering (elastic and inelastic) measurements on single 
crystals are needed to probe the fluctuations between crystal field states. El E] 
Exploring these same properties but at hydrostatic pressures below 2GPa could 
provide a valuable link to the high pressure ordered state of Tb2Ti207. [23 

In conclusion, there is strong evidence that Tb2Ti207 can not be described by 
a classical (111) Ising Hamiltonian similar to the successful model of the spin ice 
materials. 2 Spin fluctuations must be incorporated into any model of Tb2Ti207 in 
order to capture the PM correlations, and a simple description of quantum fluctuations 
via the RPA is sufficient to achieve good agreement with experiments. However, to 
study the lack of an ordering transition in Tb2Ti207, a more careful treatment of the 
quantum fluctuations is essential. Finally, we come back to the (111) Ising dipolar 
model phase diagram and note that Tb2Ti207 is very close to the AFM/spin ice 
phase boundary. jJTj There may be issues with the quantum fiuctuations, long-range 
exchange, microscopic disorder at very low temperatures, or tuning of material or 
model parameters that places Tb2Ti207 on a line or region of disordered semi-classial 
ground states. 
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